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ABSTRACT 


Measurements of temperature, wave height, and orthogonal water 
particle velocity were made in May 1973 at the NUC Tower located one 
mile off Mission Beach, California. A ee analysis of the tem- 
perature field see made using digital temperature and isotherm contour 
plots. Billow turbulence microstructure appeared to be the dominant 
mechanism with the exception of one run in which evidence of double- 
diffusion microstructure was found. Spatial correlation lengths, | 
calculated from the plot of the covariances, were of the order of 
130 cm and less when the signal was high-pass filtered for waves of 
100 seconds and longer. No depth dependence was noticed. Both fre- 
quency and wavenumber spectra were calculated and a correspondence 
between the spectra was noted at the frequency and wavenumber of the 
surface wave-induced particle displacements. 

The Thornton, Boston, Whittemore model of wave-induced tempera- 
ture fluctuations was tested and found to model the temperature spec- 
tra quite well, especially in a narrow band of frequencies associated 
with surface waves. The turbulent temperature spectrum, calculated 


as the difference between the actual and wave-induced spectra, had a 


slope near -5/3 above 0.1 Hz and more negative at lower frequencies. 
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L; INTRODUCTION 


An SHS TO RY 

For the past three years research has been carried 
out at the United States Naval Postgraduate School (NPS), 
Department of ИЛ" to study the shallow water 
interactions that occur between various oceanographic 
parameters such as temperature, salinity, surface waves, 
internal waves, and sound speed. These experiments have 
been conducted from the Naval Undersea Research and Devel- 
opment Center's (NUC) Oceanographic Research Tower located 
approximately one mile off Mission Beach, California. 
Results of the experiments are presented in NPS Theses, 
NPS technical reports and publications. This thesis is 
based on a study started by Whittemore (1973) and is an 
extension of the results obtained by Minard (1973). 

The Whittemore experiment was conducted at NUC on 8 
ann" Juüunge 19722 Measurements vere made of salinity, 
sound velocity, particle velocity, land temperature (at 
Various positions along a horizontal array) for ten pre- 
selected depths. He! concludeu that tHe relationship 
between temperature fluctuations and surface wave height 


(for swell type waves) was of the form, 


G (s) azır zan unit), 





and that a similar relationship was likely to exist be- 
tween temperature fluctuations and internal waves. He 
also found that turbulence degrades coherence between 
temperature fluctuations in sensors displaced spatially; 
however, a lack of temperature sensors prevented him 
from making definitive statements. 

Reworking the data collected during the June 1972 ex- 
periment, Thornton (1974) demonstrated a means for sepa- 
£ "ng Purbulentsand surface wave-1inducéd vertical veloc- 
ity spectral components alloving for the quantities to 
be statistically correlated. The method applied to the 
wave and velocity measurements allowed the wave energy 
density spectral components to be converted to velocity 
spectral components using linear wave theory. The com- 
puted values compared very well with the measured veloc- 
ity spectra for the moderate wave conditions encountered 
and appeared to verify the theoretical spectral transfer 
function and assumed linear system. 

Using the same data and techniques, Thornton, Boston, 
and Whittemore (1974) devised a means for separating tur- 
bulent and surface wave-induced temperature fluctuations. 
They assumed that linear wave theory could be used to de- 
scribe the surface waves; that the temperature isotherms 
were horizontal; and that the mean temperature gradient 
vasia constant. In this vay they developed a stochastic 
model describing the surface wave-induced temperature 


fluctuations. This allowed a means for calculating the 
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turbulent temperature spectrum by subtracting out the 
wave-induced contributions. The measured slope of the 
turbulent temperature spectra closely approximated the 
Kolmogorov Be q y e E a decade. ‘This Supported the 
result that the turbulent temperature fluctuations are 
statistically independent of the surface wave-induced 
velocities. 

In the present work it is intended to refine the 
analysis of the data and extend the results obtained by 
Minard (1973). Minard set out to rectify the lack of 
sensors noted in Whittemore's experiment and conducted 
another experiment at NUC on 15 and 16 May 1973. Measure- 
ments of small scale temperature fluctuations were made 
in 19 meters of water uSing a movable array of seven 
thermistors, with a short (0.015 sec) response time; and 
accurate to 0.01°C in seawater. A schematic of the tower 
and instrument mounting is given in Figure 1. In addi- 
tion to temperature, wave heights and orthogonal water 
particle velocities were measured simultaneously.  Mea- 
surements were made at various depths with either a ver- 
tical or horizontal orientation of the temperature array. 

Minard found that in the presence of internal waves 
temporal scales were’ of the order of 4 seconds and spatial 
scales were of the order of 33 cm, whereas outside the 
inrluencelof internal waves they vere of the order of 10 


seconds snd 3808 cm. It was concluded that in the presence 
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Schematic of Tower and Instrument Mounting. 


Figure 1. 
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of internal waves, correlation times were similar for 
both vertical and horizontal scales and spatial corre- 
lation lengths were similar for vertical and horizontal 


scales. 


Be OBJECTIVES 

The main objectives of this thesis with relation to 
the measurements made by Minard in May 1973 are to: 

1. Define the spatial and temporal temperature 
field. | 

p Define the spatial and temporal temperature 
Spectra and determine if an obvious relationship exists 
between them. 

3. Define the temperature gradient field and its 
variation with time to enable a comparison to be made 
between the on of microstructure with and with- 
out the passage of an internal wave. 

4. Test the stochastic model proposed by Thornton, 
Boston, and Whittemore to calculate the theoretical 


wave-induced temperature spectrum using these data. 


1.3 





LL. DHEORETICAL CONSIDERATIONS 


The creation of small scale temperature fluctuations 
in shallow water is a complex process occurring at all 
depths and scales, depending on the operative driving 
function prevailing in time and space. In the near sur- 
face regime winds over the surface produce tangential 
surface stresses both from the interfacial stress and 
from the momentum loss associated with processes such as 
breaking waves. These result in a downwards momentum 
transfer with its attendant downward flux of heat and 
salt from the surface to a depth dependent on the sta- 
bility of the water column, the gradient of shears pro- 
duced, and the wavelength of the surface waves. In this 
way the temperature and salinity in the mixed layer be- 
come virtually uniform; and, unless the stability or 


Brunt-Vaisala frequency (defined as 





N = öz , (1) 


vhere g is the acceleration of gravity, 0 is tne density 
op ! : | 

апа ҹа ls gradient of density) is less than or equal to 

zero, a thermocline develops due to the increasing con- 

trast between the properties of the water in and below 

the mixed layer. Hf there is subsitantial surface cooling, 


there may be a region below the mixed layer in which the 
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stability E is less than zero and convective mo- 
tions can develop extending to considerable depth, parti- 
cularly in polar waters (Phillips, 1969). 

The ocean bottom exerts a retarding force on any near 
bottom water moving over it. The resulting shear zone 
again causes turbulent mixing. In the absence of dynamic 
instabilities leading to a complete overturn, the bottom 
turbulence will produce a mixed layer similar to that 
formed in the surface layer and consequently intensify the 
the thermocline from the bottom. 

If a thermocline is formed, it will inhibit the verti- 
cal extent of mixing from either above or below depending 
on the depth of the thermocline and the intensity of the 
gradient within it. Turbulence produces changes which 
tend towards three-dimensional isotropy, and as such will 
be seriously inhibited by a constraint in any direction. 
In this case, the vertical motion is inhibited and re- 
sults ina flow of energy out of the other two components 
Se velocity in an efforteto boost tire vertical energy. 
Hence, Tan increase in stability as defined by an increased 
degree of stratification or density gradient may well be 
associated with the degradation of turbulence within the 
you orucolumnuğand” in particular, Hart assocľdted with 
internal waves. 

Bart debba ehe ehermocline, which tends to inhibit 


turbulent mixing from the vertical boundaries, may itself 





be a source of increased turbulence. Internal waves 
propagate along the thermocline and can lead to increased 
turbulence through shears developed between water masses 
either above or below the thermocline. In addition, 
under certain conditions the different rates of diffusion 
of heat and salt across the thermocline lead to a small 
scale perturbation in the temperature field. 

For the purposes of this thesis, the scales to which 
turbulent temperature fluctuations may be broken down 
fall into three convenient, though not necessarily inde- 
pendent, classifications. “These are: microstructure 
effects, surface wave effects, and internal wave effects. 


Each of these will be discussed separately. 


A. MICROSTRUCTURE EFFECTS 

Much evidence is now available which shows that ver- 
P1cale@maxing in they interior of Stabile fluids occurs, 
on occasion, at scales much smaller than the vertical 
extent of the water column (Woods and Wiley, 1972). 
Although turbulence can enter at any boundary of a’ fluid, 
as long as vorticity is generated, two processes have 
been proposed as capable of producing intense vertical 
mofcHositructurefacitiviity. The finst arisies from’ the diffier- 
ing rates of diffusion of heat and salt (double-diffusion 
or.ialt-£fingering). The second process is turbulent mix- 
ing driven either by surface vinds or by shear instabili- 


tives (billow turbulence). 17? occurs at the vertical 
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boundaries or on internal waves at depth. 
l. Double-Diffusion and Salt-Fingering 

This process is a local one operating where weak 
mixing has already increased the vertical gradients of 
temperature and salinity, and whose vertical scale seldom 
exceeds a few meters (Gregg, 1973). Considering only a 
two layer system involving both heat and salt gradients, 
three distinct situations can arise. 

The first is an inherently stable system where 
the upper layer is both warmer and less saline than the 
lower layer. It is of less interest because both compo- 
nents act together to produce a stable density step 
across the interface. 

The second situation is an inherently unstable 
system where the upper layer is both warmer and more 
saline than the lower layer. In a very few minutes there 
NAS a rTAáprtaRloss of heatwingthemuppcrmiayer makang it 
more dense than the lower layer, which is simultaneously 
gaining heat and consequently becoming less dense. The 
result is alternate fingers of 'heavier' water sinking 
and 'lighter' water rising, hence its name - Salt-Finger- 
MENGE 1 GUume: 2) ome Tne spattern of convection will persist 
Only to a limited depth where the descending more saline 
water Undergoes welak) overturning to form a well mixed 
layer of temperature and salinity intermediate to the 
parent waters. In such a system a series of steps will 


form beneath the interface with steps of the order of 


J 7 








Figure 2. SALT-FINGERING- The solid lines show the ini- 
tial and subsequent density profiles when the 
more rapid diffusion of heat has produced a 
potentially unstable situation at the inter- 
face. (From Gregg, 1973) 


Temperature °C 


Note: The probe is being 

| raised in the lower 
part of the record illustrat- 
ing the disturbance caused 
by the wire and the probe 
moving through the water ahead 
of the sensor. 


— Depth (m) — 





Figure 3. High resolution temperature profile of the 310-345m. 
interval in the Arctic off Ice-Island T-3 (From Denner, 
10/17: 
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20-30 centimeters and interfaces as small as two centi- 
Ne eens. 

The third situation is one in which the layers 
remain separate but where convective mixing occurs within 
each. It is characterized by a cooler, less saline upper 
layer and a warmer, more saline lower layer. Agalin there 
is a rapid diffusion of heat across the interface but nov 
it acts to increase the stability betveen the tvo layers 
vith intense convective mixing vithin each of the layers. 
In a manner similar to that for salt-fingering the convec- 
tion is self-lımiting and extends into a series of steps 
above the original interface. Neshyba, Neal and Denner 
(1971) found as many as 40 steps above a layer of varm 
saline Atlantic water that had entered the Arctic. Each 
is about 3-10 meters per step and with interfaces of the 
order of centimeters (Figure 3). 

2. Billov Turbulence 

Billov turbulence is described as free shear tur- 
bulence modified by a density gradient and initiated by 
Kelvin-Helmholtz instability (Woods and Wiley, 1972) 
3 H fe” fllüstrated in Figure 4. THis?microstructure is 
very complex in form as opposed to the rather simple, 
етти вт ılsconf 1quraittionsorschezrdensi ey steps in double- 
diffusion. 

It appears that a great majority of the flow with- 


mee ymocline is laminar and devoid of turbulent 
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Higuúre 94. BELVEN-EBEMHDOETZZSHEAR INSHABITITY 

The upper layer has a higher velocity and a lower density 
than the lower layer. The instability, generated ina 
two- Mayer, laboratory tank by rapidly raising and then 
lowering one end of the tank, caused the two layers to 
soon lose their coherence and break up into turbulent 
patches. The diagrams are based on a study by S. A. 
Thorpe of the National Institute of Oceanography in 
England. (From Gregg, 1973) 
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motions on any scale (Woods and Wiley, 1972). Occasion- 
ally, however, shear instability is triggered within 
small regions by boundary produced shears or internal 
gravity waves. In this latter case, as one of the wave 
packets travels along a sheet, shear instability is pro- 
duced at the crests of sufficient height to cause the 
local Richardson's number of the billow to fall below a 
critical value of 0.25 (Woods and Wiley, 1972). The 


Richardson number for the sheet is given by 


ghAp 


oAv” 


du x 


where "g" is gravitational acceleration, "h" is the thick- 
ness of the shear zone, "Ap/o" is the fractional density 
difference, and "AV" is the velocity difference across the 
shear. For a given value of J < 0.25 the layer will be 
unstable to a specific range of wavelengths, of which 

Lo = 7.5 x h, will be the fastest growing disturbance 
(Woods, 1968a). 

THCd.oudicion fOr termination of turbulent motions 
is that the Richardson number again increases above the 
guütiıcal value of 0.25 to about unity. When this is 
reached, the billow no longer extracts sufficient kinetic 
energy from the shear to supply the billow turbulence and 
the turbulent motions decay. Woods and Wiley (1972) found 


that in the sea the approximate scales were 75 cm for the 


Du SS ca sor thelpatch of turbulence, and 300 sec 
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Tor tremwlifetimemwot the sturbulénce, 


D SURFACE WAVE EFFECTS 

During the past few years considerable effort has 
been put into understanding the effects of wave motion on 
the thermocline but relatively little has been done to 
understand the internal temperature effects such waves 
produces Clearly, temperature fluctuations in the near 
surface region can be caused by a number of factors, in- 
cluding surface waves and their associated turbulence, 
internal waves, and advection. Surface wave effects on 
temperature fluctuations are usually considered only to 
depths of less than one-half of one wavelength due to 
the exponential decay of particle motion with depth. How- 
ever, the energy densities of surface wave-induced fluc- 
tuations are in the same frequency band and of the same 
order of magnitude when they are present simultaneously 
and hence are very difficult to separate. Therefore the 
spectral separation of the turbulent and wave-induced 
temperature fluctuations will be considered in greater 
(ea ll later in this thesis. 

The temperature field can be considered as a turbulent 
fadla navıamng som "icnergy\ containing negion, fan "inertial" 
T and an "energy" dissipation region analogous to 
that hypothesized for velocity fields. In addition, since 
the nates of thermal and viscous dissipation are different 
in the ocean, the smaller scale temperature field will 


Bends te persist longer than the velocity field. 
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E süz b Ca ve ndüceki Temperatume EluUctuations 


Earlier, mention was made of a paper (Thornton et 
al, 1974) setting out a relationship between wave- induced 
temperature fluctuations and the waves. The major portion 
of what follows is taken directiy from that paper and is 
included here, in a somewhat abbreviated form, for com- 
pleteness. 

Thornton (1974) had previousiy shown that linear 
wave theory estimated the wave-induced water particle 
motion very well under small amplitude wave conditions. 
Since it worked for particle motion, linear wave theory 
was assumed to estimate the wave-induced temperature fluc- 
tuations. This implies that the flow is irrotational. 

It was further assumed that the temperature fluctuations 
are small, allowing density to be considered constant over 
short distances. Hence, temperature is considered to be 

a passive scalar with buoyancy effects neglected. Assum- 
ing temperature is a conservative property, at least in 
bilie @bodyvtor.thé.flow”the conservation of temperature 


flux Iis given by 


90 90 EL 
— + u. <— + v — = 2 
t “ox öz = E 
where O(x.,z,t) is the measured temperature, u (x.,z,t) 
are the horizontal velocities, v(x, ,2,t) mothe vertical 


VeELGETUY. x: arc Eu" horizontalscarues2uoanseoordinates, 


z 1520006 060 0731 cartesian coordinate: The equation 
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considers the individual time change, advection, and dif- 


fusion of the passive scalar temperature. In terms of 
their means, turbulent fluctuations and wave-induced 


fluctuations the quantities are defined as 


f 
Il 
DI 
$ 
о 
4 
о 


where the overbar indicates a time-averaged mean, the sub- 


script w indicates contributions from surface waves, and 
the prime indicates turbulent motions. 
Assuming that 
(1) turbulent and wave-induced motions are sta- 
tistically independent, 
(2) mean temperature isotherms are horizontal, 
(3) mean velocities are zero, and 
(4) the mean vertical temperature gradient is 
constant, 
allows separation of the turbulent motions and Equation 


(2) reduces to 


20 20 2 x 
Da ә ив 2 Di 


where the velocities and the assumption of irrotational 
flow are as specified in linear wave theory. 
Temperature fluctuations caused by wave-induced 


5500100 16 кс a priori assumed equal to zèro and 
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therefore 


90 90 


(n) i IX. È dər 2 5” 





This leaves a first order, constant coefficient, 
homogeneous differential equation, if the mean tempera- 


ture qradient is constant, and therefore 


. (4) 


This is easily integrated to give 


98 sinh k(d+z) 


77: E a 


İG cos (kix;-0,t) (5) 
where V (ži Zt) is the wave-induced vertical water-par- 
ticle velocity specified by linear wave theory, 'a' is 
the surface wave amplitude, k is the wave number, G is 
the radian frequency, d is the depth, and z is the nega- 
tive downwards displacement from the mean water surface. 


(0) 


im terms) of=the®@sipectral™transiier function, Ho 


ehe kiuare bracket in Equation (5));, and the#sinusoldal 


surface displacements, n(xi,t), Equation (5) reduces to 
О (x, z, t) = Hg (0) n(xi,t) (6) 


Equation (6) satisfies Equation (3) exactly, which means 
that moderate swell type waves do little or no mixing, 


Dursonivzsserve. to pump the thermocline Up and down. 
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2: Turbulent POSE ic IG 


In a manner similar to that discussed in Thornton, 
MUST EN ema Der bv Thornton, Boston, nd Whittemore (1974). 
discusses the separation of the turbulent temperature 
spectrum from the wave-induced temperature spectrum. For 
completeness the arguments are repeated here, 


The wave-induced temperature spectrum, So (О), ls 
w 


calculated from the wave spectrum, nr using the com- 
plex transfer function, Hg (0) (given inğEqguatıon (6)); in 


the form 


2 
o (О) = [Hg (0) | 06022 (7) 


S 
” n 


In terms of the total temperature spectrum the 


turbulent temperature spectrum is given by 


Sao es s D api əə : (8) 


If the cross-spectrum is expressed in terms of 


its co- and quadrature spectra 


EE 100 (0) 
n n n 


then Equation (8) becomes 


Sg. (0) 7 Sg(90) + ə - 2H, (0) su (9) 


which gives the turbulent temperature spectrum in terms 


of measurable quantities. 
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Assuming statistical independence between turbu- 
lent and wave-induced temperature leads to further sim- 


plification, since 


(G) = Sg (0) 


and 


2 
E ER O E (G) = | Hg (0) | S, (0). 
vn v 
This allovs the turbulent temperature spectrum to be 
given simply as the actual temperature spectrum minus the 


vave-induced temperature spectrum, or, 


2 
Sg. (0) 7 SQ(6) - [|Bg(o) | 5, (0) . (10) 
Further the coherence between temperature fluctua- 


tions and the waves is then given by 


SOJA 
Te 7-5” T D 
n 0 


Cz INTERNAL WAVE EFFECTS 

Internal waves are subsurface waves existing between 
layers of different density or within layers where a ver- 
Fıcalımdensiıty gradıcnc sexists. "They occur often due to 
a CEN CARO causes such as flow over an irregular bottom, 


atmospmere disturbances, tidal forces, and shear flow. 
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The theory of the existence of internal waves implies 
perfect coherence of plane waves Over space; whereas in 
practice naturally occurring internal waves are not so 
coherent. In addition, the ocean contains a number of 
generating sources which may reinforce or cancel each 
other (LaFond, 1966). 

In shallow water, at the same site as this experiment, 
LaFond (1966) found internal waves of amplitudes ranging 
from 2 to 238felet®with mean*amplitude*of 5.6 feet. Win 
general, the magnitude of the shorter period waves was 
found to be inversely proportional to the gradients in 
which they were found. Theoretically, internal waves 
have periods ranging between the lower inertial and higher 
Məehisala or istability frequencies. LaFond found periods 
ranging from 4 to 10 minutes superimposed on longer diurnal 
cycles and four and one-half day cycles. He also found 
that internal waves moving towards the shore had an average 
speed of 0.31 knottsMand ranged from 0,11 to 0.60 knots. 

Physical intuition would lead us to suspect that tem- 
perature fluctuations associated with internal waves should 
exhibit a linear relationship between the temperature 
gradient and the wave amplitude. Therefore a reasonable 
expression for the temperature fluctuations associated with 


internal waves would be 








where the overbar indicates the time-averaged mean and 
A (t) is the instantaneous amplitude of the internal wave. 
The temperature fluctuations due to internal waves 

are usually easily identified because of the narrow range 
of frequencies where significant energy is present, and 
the presence of a valley, in the spectral energy density, 
separating the surface and internal wave regions. Further, 
energy contained in the turbulent temperature fluctuations 
is small relative to the energy of the internal waves, and 
1s, therefore, only a minor feature at these low frequen- 


cies. 
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LUL: EXPERIMENT AND DATA ANALYSIS 


A. EXPERIMENT 

Measurements were made at various depths with either 
a Vertical or horizontal orientation of the temperature 
array as summarized in Table I. 


TABLE I. Summary of Experimental Set-up for Minard's 
May 1973 Experiment 


Depth to š 
Run top of Array Internal 
Number] Date frame(m) | disposition waves 


vertical 
vertical 
vertical 
vertical 
horizontal 
horizontal 
horizontal 


horizontal 





id 


Faqure#5 iliustratels the array configuration for the 
first four runs, conducted on 15 May, where the placement 
BrgtBesthermistors was vertical. igure 6 illustrates 
Eheonmcontiquratiıon of thegarray jon the second day, 16 May, 
when the thermistors were in a horizontal line. Rune 
was conducted in T configuration but in all subsequent 
runs the array was rearranged such that Thermistor 1 was 


Tom 777777 1 accaumc4 inches (60.8 cm) below Thermistor 5. 
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Figure 5. coni igon of Array for Runs 1, 2, 3, 4. 


BL 













o 24.2 cm 
\ | Wave Guage 








te 
O 
di 
o 
€ 
Š 
no 21.5сп 
g lees | 
AN © < N "О N 
Wu. Uno "ës = == | == == 
2-1: co e 
NO o o A 
— e © “О 
po 
Thermi 
J seaward 
Elquroao lltonfıigquratlon of Array for Run 6. For Runs 


7 and Ba, Thermistor 1 was removed and placed 
60, Seemebelow  Whermistor 5. 


32 








This was done to give an indication of the vertical tem- 
perature gradient. Data were recorded on a Sangamo 3562 
analog tape recorder using F.M. electronics. These data 
were played back at 60 ips (32 times speed-up) and digi- 
tized at NPS on a hybrid (analog-digital) computer, CI5000 
EXDS 9300. The sampling rate of the data analyzed was 
reduced to 6.25 Hz leading to a Nyquist frequency of 3.1 
Hz. This reduction was not considered critical, since 
the spectra became essentially flat and the signal lost 
in noise before the 3.1 Hz frequency was reached. Data 
were analyzed on the school's IBM 360-67 computer. The 
environmental conditions consisted of light winds, small 
amplitude waves and swell, low overcast skies, and occa- 
sional internal waves during the two day period. Further 
details of both the experimental set-up and the analysis 


Rarticulars can be found in Minard (1973). 


B. DATA PREPARATION 

Initially all the signals were graphically displayed 
on an eight channel strip recorder for visual inspection 
ufutncirecorde, pofTarlity/Terror and anomaly checks. Con- 
firmation vas made that runs 5 and 8 vere unusable for 
digital purposes. Run 9 vas designated Run 8a and re- 
placederuns 79:07 analysıs Durposes” since they vere” made 
at the same depth vith the same configuration and in 
nearly the same time frame. The analysis was carried out 
using the nine-track digital tapes previously prepared 


BEOM EMG analog records by Minard (1973). 
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C. METHOD OF ANALYSIS 

The variances and covariances were calculated from 
the time series of each run by first removing any long 
period trend and any D.C. bias component. In addition, 
some of the time series were digitally high-pass fil- 
tered to filter out waves of period greater than 100 
seconds prior to computation. Spatial correlations were 
computed by averaging the covariance (C.V.) between the 


first thermistor and each of the others -such that: 
CV — Jm, tex) Tn, (e ,x+0x,) at 


where T is the length of record (= 20 min). The C.V. 
was then plotted against the appropriate separation dis- 
tance along the array using the first thermistor as a 
reference. The variances and covariances were normal- 
ized with respect to the variance of the first thermis- 
Bor in thewarray. Both the normalized and the unnormal- 
ized spatial correlation functions vere Fourier 
transformed to obtain the normalized and unnormalızed 
wavenumber Spectra for the runs. 

Next; the gradients present in Runs l, 4, 7, and 8a 
were calculated since they were to be used in the theo- 
retical wave-induced temperature spectrum calculation. 
First, the digitized signals had to be re-converted to 
erue “Cnperature” since during the recording and digi- 


tizing processes a number of gains and biases had been 
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introduced. The original temperature recording had been 
a voltage conversion of the temperature variation around 
a balance resistance in the wheatstone bridge circuit of 
each thermistor. Therefore, the actual temperature 


would have been of the form 
T= T + AT , 


where the overbar indicates the mean temperature and AT 
is the variation about the mean. 

In San Diego, all signals were amplified before re- 
sording: Therefore the recorded signal would be G1*AT, 
where Gl was the appropriate gain. Еплкет пи Һе ал јл 21 21 ПҸ 
process a bias (В1) was applied prior to the computer 
gain (G2), which was wired into the analog circuitry. The 
digitized signal, in terms of the original fluctuation, 


would be 
D.S. =- G2(G1. * AT +-Bl) . 


The original temperature in terms of subsequent biases 


and gains would be 


– ( D.S. - Bl”G2) 
= T + , ——— —— . 
e Ge s. G2*Gl 


The gradients were then calculated for the selected 
runs using the temperature difference between thermistors 
il onun, vhichnad o constant Spatial separation of 24 


mehes O Gm), as representative of the whole array. 





A separate program was written to perform the spec- 
tral analysis of the time series. The data were read 
into files corresponding to Thermistor 5 and the Baylor 
Wave Gauge. The records were then detrended and the 
mean removed. Autocorrelation functions for each file 
and a crosscorrelation function between them were calcu- 
lated. The autocorrelations were normalized with res- 
pect to their variance and the crosscorrelations with 
respect to the square-root of the product of the indi- 
vidual file's variances. From these functions, the auto- 
and cross-spectra were calculated by Fourier transforming 
the respective temporal correlation functions, after 
application of a Parzen window filter to account for the 
finite length of the time-series. 

From the Baylor Wave Gauge spectrum, the wave-induced 
temperature spectrum was calculated uSing the spectral 
Budusner irunctiOnuderived"in Egquation” (5) In addition, 
coherence and phase angle were calculated between the 
temperature and waves. The theoretical turbulence spec- 
trum was calculated from the difference between the 
actual temperature spectrum and the wave-induced tempera- 


Eure?”speçcçtruüm as güven in Hquatıon (9). 
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IV. RESULTS 


A. ANALOG TEMPERATURE RECORDS 

Examination of the analog records of temperature, re- 
produced in Figures 7 (Run 6) and 8 (Run 1), show that an 
internal wave was present in Figure 7 and not in Figure 8. 
The common features of each are the high levels of tempera- 
ture fluctuation and the general consistency of prominent 


features throughout the temperature array. 


B. TEMPERATURE. FIELDS 

Some long term fluctuations were noticed on all runs 
and four runs (2,3,6,7) show significant internal wavie 
aetivity. 

| Runs l and 4 show patches of increased turbulence 

which may or may not have been remnants of internal waves. 
In any case, the scale of the turbulence was such that 
the whole of the vertical array was affected. Figure 9 
пар üZuUngğtne digitized data of the first 15 minutes 
of Run l with a data point plotted every 4 seconds. The 
total vertical extent of the array covered 1.625 meters 
and over such a short distance the temperature difference 
between thermistors was never greater than 0.5°C and 
usually less. 

Figures lO, er pilot or the isotherms for Run l, vas made 


DHe onr ert ing each thermistor signal to an actual 
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Figure 7. Analog Temperature Record-Run 6. 
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Figure B. Analog Temperature Record-Run 1. 
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temperature reading in the manner described in Section 
III-B. Subroutine "CONTUR", from the Schools IBM/360 
source library, was used to plot the isotherms at inter- 
vals of 0.3°C with a data point every 4 seconds of real 
time. The temperature data from the time series was 
linearly interpolated every 3 inches across the array, 
at each time increment, aS input to the subroutine. The 
plot shows a relatively warm layer for the first 15 cm, 
underlain by a relatively cold layer in the next 15 cm, 
which in turn is underlain by a relatively warm layer 
covering the rest of the array. 

Although no measurements were made of salinity in the 
detail required to solve this question, the salinity 
gradient must be slightly positive or at least isohaline 
bor the first 30 cm df the' array, tand then increasingly 
positive to at least 60 cm in order for this system to be 
dnamıcalıy Seable, L£ this)is ehE€ case, then the possi- 
511:677407 double—difTusion type microstructure, with 
SAtablleMstratitication fand convective' mixing within each 
layer, exists between the central and lower layer. The 
increased gradient of temperature between 15 and 30 cm 
is consistent with an inherently stable double-diffusion 
mechanism. TA feci the salinity gradient is negative 
ити Һә ис ws Omcm then the possibility of salt-flingering 
exists (inherently unstable double-diffusion), although 
smeresresno direct evidence ın the 1sothermifield to sug- 


ОР Си mechanism rs occurring. 
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The decrease in the temperature gradient in the lower 
half of the array (Figure 11) may be an artifact produced 
by the array spacing and the subsequent interpolation of 
Polines in them ontour plot. 

The blobular formations within the warm upper layer 
and in the cold central layer are probably due to inter- 
actions with small amplitude internal waves which were not 
recognized as such from the temperature records alone. It 
would appear that the layers are being alternately 
squeezed and relaxed with a periodicity approximating 
3 minutes. These blobs compare well in time to the in- 
creases and decreases in the turbulent temperature evident 
in Figuře"9, 

Runs which show significant internal wave activity are 
itemized in Table II by run number and "periodicity" 


measured from the analog records. Periodicity is at best 


TABLE II. Internal Wave-Periodicity i ` 


run number period (sec) 





approximate due to the limited number of internal waves 
pesi Conai.) quegli lisa plot of the first 15 minutes 


of Run 3 showing the temperature variation for Thermistor 1 
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5777 7 emar ray mate 6010 - meters belov the Surface) and 
Thermistor 6 (tovards the bottom of CS array). One in- 
ternal wave had just finished at the start of the plot. 
A second was starting at the end of the plot and con- 
7 nucceefuto Kigüreğl2, thewmnexe 15 minutes of Run 3. Tt 
illustrates the periodicity of the internal wave. Figure 
13 is the computer plot of the isotherms for the first 15 
minutes of Run 3. The intrusion of colder water into the 
upper portions of the array and its subsequent withdrawal 
as the internal wave passes is quite noticeable. The 
salinity profile corresponding is available but is 
probably a positive or isohaline gradient with increasing 
depth. Double-diffusion type microstructure would not be 
present under these conditions of stratification and 
billow turbulence microstructure would be more likely. 
The thermistor array during Run 3 was located just 
below a slight thermocline which had been decreasing in 
depth and, to some extent, intensity from the time of 
Run l(Minard 1973-Figure 1). This suggests that the 
presence of internal waves in the period of time from 
Run l1 to Run 3 was contributing to the depression and 
erosion of the thermocline. Figure 13 shows the extreme 
variability and "fluctuation in an individual isothermi. 
Somes Ol@theslluctuation 1s no doubt attributable to tur- 
bufcnceeğuvhicn could Defa cause of thai thermocline erosion. 
However, the assumption that the fluctuations are solely 


caused by turbulence must be tempered to some degree 
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since successive peaks in the fluctuations are separated 
by about 12 seconds, a figure that is suspiciously close 


to the surface wave period. 


cs SPATIAL CORRELATIONS 

Table III contains a listing of the covariances be- 
tween Thermistor 1 and the others, and the individual 
variances Of teach thermistor, for all of thie runs for both 
non-filtered and high-pass filtered data. To accomplish 
the high-pass filtering, a Fourier analysis was made of 
each 20 minute segment of data; the Fourier coefficients 
generated were used to generate a time series describing 
mies waves of 100 seconds and greater. This series was 
then subtracted from the original time series which was 
then the filtered time series used for analysis. 

One of the most common ways to define a scale length 
of some inhomogeneity is to measure the length required 
forthe correlation to fall to Wir (0.368) off TES original 
value. The correlation functions, shown in Figures 14 and 
iS = cre of a very long scallesteature and only one run in 
eRchirigure’ actually falls to 0.368. Thismphenomena is 
inteorpreteddlas being due tö the extremely long scalle#of 
the internal waves and their attendant higher overall 
energy. Consedgquentlysrtowas decided to filter out the 
waves er periods 100 seconds and longer. This was accom- 
plished daigıiıtaliy, as lexplainedi/flarlier. The results, 


meer hem tal £[ilvering of the temperature signal, 
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are illustrated in Figures 16 and 17 where it is appar- 
ent that waves of 100 sec and longer had been contribut- 
ing substantial energy as evidenced by the overall 
lowering of the spatial correlation curves when filtered. 
Table IV lists the correlation lengths for each run 


after the original signal had been high-pass filtered. 


TABLE IV. Correlation Lengths 





Run Depth-Surface to Correlation 
Number Top of Array (m) | Length (cm) 
1 
2 
3 
4 
6 2.7 40.0 
6.2 | 50% 0 
8a Iı l 21.4 


Although it might be expected that the shortest cor- 
relation lengths would occur at the surface, where waves 
dominate the structure, no such relationship is apparent. 
Herhapsi thisf can be explained by the fact that” thelarray 
dimensions are small compared to the depth and breadth of 
enhe surrounding kocean arĉa; Or Enetscale of the turbulence 
may be quite large compared to the depth of the water; or 
maybesıt further aldllustrates that waves do little mixing 
but omly pump thie walter up and down. 

BeconsäistentestrendZwäisenöticed in that the horizontal 


drray always had shorter correlation lengths than the 
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vertical array at similar depths. However, this result 
is noe conclusive since horizontal and vertical array 
runs were made on different days. 

Further examination of Table III shows that the runs 
with significant internal wave activity have higher 
energy as evidenced by the higher numerical values of 
variance and covariance. This trend is present even 
when the signal was filtered for waves of 100 sec and 
greater except that Runs 2 and 3 are not significantly 
greater in energy than Runs l, 4, and 8a. This indicates 
that Runs 6 and 7 contain significantly more energy at 
all frequencies than Runs 2 and 3. Also the difference 
in covariances between signals not filtered and filtered 
show that Runs l and 8a contained much less energy at 


periods greater than 100 sec than any of the other runs. 


D. WAVENUMBER SPECTRA 

The normalized wavenumber spectra were calculated by 
Fourier transforming the spatial correlation functions 
(Figures 16 and 17) calculated as explained in Section 
len. Be Thesspatial correlation was calculated “at equal 
inerements@aliong the’ llength of thel array by Minearly in- 
PBerpolating between the calculated points on the array. 
Figure 18 is a composite of the normalized wavenumber 
spectra when the array was vertical (15 May); and Figure 


NOP 15) ağcomposite when the array waisi horizontal (16 May). 
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Figure 19. Horizontal Array Wavenumber Spectra. 
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The spectra show a high degree of consistency between 
runs. In the vertical array, Runs 2 and 3 show slightly 
more energy at all wavenumbers than Run 4 but less energy 
tuangSRuns ls It would appear that Run 4 had a major por- 
tion of its energy at periods greater than 100 sec as was 
the case for Runs 2 and 3 with the internal waves. In 
the horizontal array spectra, all runs are very Similar 
in shape whereas the vertical array Spectra varied con- 
Siderably in shape. It is felt that this may be inter- 
preted as a difference in the source of turbulence on the 
two days. 

Each of the spectra has what is believed to be an 
anomalous spectral energy increase at higher wavenumbers. 
Gregg, Cox and Hacker (1973) observed a similar increase 
at slightly higher wavenumbers in the temperature gradi- 
ent spectra. They attributed the initial rapid fall-off, 
below 0.06 radians per meter (rpm), to dynamical processes 
dominated by the nearly exponential mean temperature 
gradient with depth. For wavenumbers greater than 0.06 
buerziesisztHhanz3T.47apmäthey felt the spectrum was domi- 
na aed by thoas inregulanliy spaced high gradient regions. In 
the reqionmef@interest, above sal.4 rpm, where high levels 
Or microstrúücture exist the spectrum: takes the form of a 
minimum followed by a maximum at higher wavenumbers. This 
last region is in general the same wavenumber space as the 


increase in our spectrum. Whether the increase is due to 





the same cause is uncertain; but there were high levels 

of microstructure activity in all the runs. However, it 
may be that this increase is due to other factors such 

as the manner in which the correlation function was calcu- 
lated in deriving the spectrum. No attempt was made to 
fit the correlation function to other interpolation forms, 
e.g., Gaussian or exponential, near the origin since the 
structure function when calculated did not show that the 
time series were particularly stationary. 

It seems reasonable that the wave-induced contribu- 
tion to the wavenumber spectra would occur over a narrow 
frequency band corresponding to the frequencies of high 
coherence in the frequency spectra, Consequently, the 
Spatial scale in wavenumber space should show corres- 
pondence to the orbital excursion of the wave-induced 
Particle motion within this narrow band. The scale of 
the particle displacement is calculated using the root- 
mean-dquare (rms) wave height in each run. The vertical 
displacement, &(t), using linear wave theory, is related 
bomcherwaverproßıle, p (t), Dy 


shnh k (zd) 


-. . Lex kd 


lt) 


where the term in brackets is the spectral response func-' 
tion, Н. (О). The maximum displacement can be approximated 
at the peak in the frequency spectrum as 


52 : rms 
E ELLA Ls (12) 


90/92 
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where H (€) is the response function at the frequency 


max 
of maximum coherence and ga is the rms wave height. 


The horizontal excursions can be similarly calculated, 


using linear wave theory, as 


IS [E ax] coth k(zta). 


The results of these calculations are presented in Table 







V. 
TABLE V. Particle Displacements 
Wave- Wave- | Expected 
Depth |Freq. | & num- с num- wave- 
. ber əə ber number 
Run (Hz) (m) (rad/m) (m) (rad/m) (rad/m) 


0 3.6 0:092 120.76 Вә 12253, 5.0 EZ 
0298 6.4 E. OSO 8.0 
0.66 Ss 2520017750 8.0 


8a 1.286 13.8 0.092140. 30 821.1 1.41] 4.5 1220 


These scales are representative of the diameter of the 
nahber”"uartrıcle orbıt. It would ‘be expected that 86 there 
"soda neak in the frequency spectrum (Figure 24) corres- 
ponding to surface wave-induced motion, then there should 
be a DeL in the wavenumber spectrum assoclated "ICH tHe 
same motion at the wavenumber of the excursion. Examina- 


tion of Figures 18 and 19 shows that although there is no 
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peak in the spectra there is a definite change in the 
slope of each spectrum in the vicinity of the scale 
length" calculated throughs Equaltion (18). The fact that 
the unÜersections do not lie exactly at the change of 
slope may be due to the use of the rms wave height 
rather than the median height or some other value. Also 
as would be expected, the scale length, £, decreases as 
the depth increases illustrating the decrease of the 


orbital motion with depth. 


E. GRADIENT FIELDS 

Gradients for Runs 1, 4, 7, 8a were calculated and 
Sheydiqrtal data plotted. Table VI isal compilation of 
the calculated mean gradient, the coherent gradient used 
in calculating the wave-induced spectrum, and the gradi- 
ents as determined by Minard from the bathythermographs 


IMrubmd, 19751. 


TABLE VI. Gradient Data 


Top of the Top of the BETI Calculated | Cohereñt 
Run thermo- frame gradient gradient gradient 
Number {cline (m) (m) 57m) ( "C Am) (C m) 


5 2009 
0.19 notscale. 
not calc. .096 


-0.41 -.07 


0.10 





02 








As can be seen from Table VI, the gradients taken 
from the Bathythermograph (BT) are in general agreement 
with the time-averaged mean gradient from the time-series, 
with the exception of Run 7. The difference can be 
accounted for because of the difference in length over 
which the gradients were determined. The gradient deter- 
mined from the BT's was originally averaged over a dis- 


` 


tance of 18 meters and was all that was available, whereas 
the calculated mean gradient was determined "m 24 inches 
(0.6 m), the distance between Thermistors 1 and 5 on all 
runs. The difference in length over amen the gradient is 
calculated is important and can be illustrated from the 
results of Run 1. The mean calculated gradient was deter- 
mined to be 0.93°C/m (Table VI). Figure 9 shows the tem- 
perature profiles of Thermistors 1 and 6 of Run 1 (same 
time frame but Thermistors 1 and 6 are separated by 14 
inches (0.30 m) more than Thermistors 1 and 5). Here the 
maximum temperature difference is about 0.5°C, and since 
the separation was 38 inches (0.96 m) between thermistors, 
the gradient would’ be atmaximum of 0.523°C/nm. It clan bie 
seen from this that a small change in the length can re- 
sült inga largefchangeg in the gradient. 

Mü canğafso biel. seen that thie ¿BT amd calculated gradi- 
ents?are much greater than the gradient required to match 
coherences in the wave-induced spectrum. This indicates 
that the spectral model would overestimate the wave-induced 


temperature fluctuations if the actual mean gradient were 


6 3 


Lo » - 





used. Thornton et al (1974) noticed the same result 

and calculated the effect of buoyancy to try and account 
for some of the difference. They concluded that "buoy- 
ancy effects do not dampén the motion, but serve to 
amplify it relative to the homogeneous case." Hence 

they could not explain why the mean gradients are greater 
than the coherent gradients since the wave solution would 
be improved by less than 1% by including the buoyancy 
effects. What is obvious is that the specification of 

a constant mean gradient is not completely realistic in 
the presence of internal waves and turbulence. 

Figure 20 is a plot of the digitized data of Thermis- 
ors Trandzs during the’first 22 minutes of Run 1. 
Thermistor 1 was 0.6 m above Thermistor 5 and therefore 
al posa tivefgradiemt Ws indicated, l.€., temperature de- 
crease with depth. Figure 21 shows both the gradient and 
Thermistor number ] temperature plotted together for the 
Bırsrs22smınutesgor Run*l¿ Each point On the plot 15 the 
averaged gradient over 1.44 seconds. Thermistor 1 is 
necilotrtted"iromsFıqoüre 2Oebuteatmwam@slightly different 
scale. The arrows labelled A kind B in both Figures 20 and 
21 "serve)b to indicate noticeable points of gradient incréase. 

Figures 22 and 23 show the same information, as ex- 
pHarned above, tor §Run 7 vhene? Thermistor 5 is now 0.6 m 
abboveğThermustor ls Again a point is’ marked (C) to illus- 
Bratesansıvncreasesinegräädient> What is interesting, how- 
ever, is the difference in the temperature profiles and 


gradients between Runs 1 and 7. 
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Figure 21. 
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In Run 1, the temperature profile appears to be as 
'turbulent' and perhaps more so than Run 7. Run 7, how- 
ever, is marked at the beginning by the passage of an 
internal wave. The gradients in Run 1l are relatively 
constant, being uniformly positive with a deviation of 
less than 0.5°C/m from the mean. On the other hand, the 
gradients in Run 7 are extremely variable, being both 
positive (indicating stable conditions) and negative 
(indicating unstable conditions) with extremes ranging 
from +0.58 to -1.04°C/m. In addition, it appears that 
most of this record is best characterized as dynamically 
unstable. This is certainly indicative of intense turbu- 
lence associated with the internal wave and there may be 
some indication that the wave has 'broken'. In the area 
Se point D on Bligure 23, the gradient is almost zero for 
a period of nearly 150 sec and represents a zone of 
nearly homogeneous water brought about from the intense 
turbulent mixing. In an analog to surface waves, this 
zone could be termed an 'internal swash zone'. The calcu- 
Wc ed gradients for Rumse4 land Bagareesimilareto the re- 


sülüs for Run l. 


5. ZTEMPDEPALURL EFELOUENOECYT SEM RƏ” ANALYSLS 
mu Us of Rune o Wa. Li be discussedgmore fully 
Elan chegotherzrunsSsiinger it most cHearly demonstrates 


the effect of surface wave-induced fluctuations. The 


69 





other three runs analyzed (1, 7, 8a) will be discussed 
only insofar as they differ from Run 6. 

T Surface Wave-Induced emporcturocdriuctuations 

Equation (7) vas used to calculate the vave-in- 

duced temperature spectra. Initially the mean tempera- 
ture gradient, as calculated from the time-series, vas 
used in the transfer function given in Equation (5). 
This led to an overestimation in the wave-induced tem- 
perature spectra. The variability in the temperature 
gradient A lavas considered to Bömainliyan function ofi 
the internal waves present and consequently the use of a 
mean gradient was hardly realistic. If an approximation 
to a mean gradient had to be used, it may have been 
better to use the median, or most commonly occurring, 
gradients Since a statistical analysis of the gradient 
was not possible it was decided to use the method of 
establishing a coherent gradient as Thornton et al (1974) 
had done. Therefore, a gradient was used which matched 
the measured peak coherence value between waves and tem- 
perntumse withga calculated (theoretical) coherence using 
the calculated wave-induced temperature spectrum substi- 
Cu rtedismtoMEdguationstil)- their rational for doing this 
was that if the turbulent and wave-induced temperature 
fluctuations are statistically independent then the 
measured coherence should be closely approximated by this 


Sod tomer rovldGadethestranscer function was correct. 





Since the model is based on linear theory the only "free 
variable" in the transfer function is the gradient. This 
coherent gradient (Table VI) is Significantly different 
from the mean gradient. Table VII shows the appropriate 
values of depth, frequency, coherence, temperature spec- 
tral density, wave-induced spectral density, and average 
phase angle (between the wave and temperature spectra) 

at the point of maximum coherence in the wave and tempera- 


EE EC ms pe dure. 


TABLE VII. Frequency Spectra Data 


Temp. W-I Temp. | Average 
Coherence | Spectral Spectral Phase 
max. Value Value Difference 





Figure 24 shows the measured temperature and calculated 
wave-induced temperature spectra, and the coherence and 
phase angle between waves and temperature which are dis- 
eussed in the! following sections. 
Zn CO Ne Hen cle 
The 1555: coherence between waves and tempera- 


ture is shown in the middle graph of Figure 24 as the 


solid Tife. The (theoretical) coherence between waves 
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and wave-induced temperature was calculated from Equation 
(11) assuming statistical independence and is illustrated 
as the dashed line on the same figure. It is only in the 
relatively narrow bandifrońňm-0.05 to 0.3 Hz that there is 
significant coherence between the spectra of waves and 
temperature. This corresponds to periods of between 3.3 
and 20 seconds which is most certainly characteristic of 
surface waves and swell. The theoretical coherence, 
which was adjusted to match the measured coherence at the 
first peak value, approximates the measured coherence 
very closely. Hence the separation of the two spectra 
can be interpreted, from the increasing lack of coherence 
with increasing frequency, as an increasingly high ratio 
omt turbulence (noise) to coherent surface wave-induced 
temperature fluctuations. 
3. Phase Difference 
The last graph in Figure 24 is the measured spec- 
tral phase difference between the waves and the tempera- 
ture. A visual inspection of the figure reveals that 
over the range of significant coherency, the average phase 
“ru Horoncob ic” 153: 797 c1ues Toc the Pother 9cuns areğlısted 
in Table VII). In theory, the phase of the waves should 
lead the temperature by 180°. The difference was un- 
eXrpectedseincesthezsresults of Thornton et al (1974) indi- 
cated a much closer relationship between theory and actual 


measurement. Lu Suda Hey found that the Span over 
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which the phase difference was close to 180° was much 
wider than the frequency band of high coherence. This 
occurred in these results only once, in Run 7. Also, 
the phase difference appeared to be depth dependent be- 
cause in Runs 6 to 8a where the depth of the sensor 
(Thermistor 5) increased, the phase difference decreased 
I Bon, 1 55 to nl 059, 

The reasons for these differences are not readily 
Tontheoming. Mm Run 1, as: noted: earlier, the stratifica- 
tion was different from the other runs in that over the 
length of the array the temperature field was highly sug- 
gestive of some double-diffusion mechanism operating. 
This may suggest that in such cases buoyancy plays a 
much greater part in the destruction of coherence between 
waves and temperature than in cases where temperature 
steadily decreases with depth. 

The actual difference between the average phase 
Ci Hference kand 1802 wise even more puzizling. With their 
sensor at 5.5 m Thornton et al (1974) found almost exact- 
Nc LL. 5.0 O Ts ela £ fe Hences Ecom about 2,05 to 0.22 Hz. 

The results of this experiment at a lesser depth indicate 
SMphesefchangefetf about 155°" (about 25° dijifierence) which 
decreased as the sensor depth increased. At first the 
difference was attributed to the horizontal displacement 
daz 5 5“75oho2s1ag 7f3Thermistorc 5 "from the Baylor 


Wave Gauge. The first order expressions for the wave and 
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temperature fluctuations are: 


a cos(kx - Ot) 


= 
il 


од H (0) a cos (kx - Of) 


= “H, (8) a cas [180° = (kx = OE)] X 


Using the wave gauge as a reference position, the phase 


angles become, 


2 = -0t 
n 
Le = 180° - (KAx - ot) , 
ScHkehatsthe plläse difference, A2 = 2a 2 , is 180” #£Mx 


0 n 


is zero. The horizontal displacement was 0.19 m which 

TOL thws) calculation jat the frequency of maximum coherence 
qclded?a differenc€ of 0.52° from#180°®. This Ys a 
derınlitessndicatıon that "all is not as simple As! the 
assumptions in the model would have us believe, En cananı 
waves, which of themselves would be 180% out of phase with 
che temperature fluctuations, would certainly complicate 
the picture because of their interaction with the surface 
waves and the vertical displacement of the thermistor and 
the wave gauge. Buovaney Grrects Erougqne”ubourdoy small 
Scale instabilities such as those noted in Run 7 may be 
AMO Cer mictor mend Lever tire’ Gause, ЈС сап be Said” from 
emesemresultes that the phase dirferenee#is not due to the 


surface wave-induced motions alone. 





4% Turbulent Temperature Fluctuations 


If in the consideration of the theory we assume 
that turbulent and wave-induced temperature fluctuations 
are statistically independent, then the area between the 
spectral energy density curves on Figure 24 should repre- 
sent the turbulent energy as given by Equation (10). 
Figure 25 is a plot of the turbulent temperature spectrum 
for Run 6. Thornton et al (1974) calculated the turbulent 
temperature spectra considering the possibilities that 
turbulent and wave-induced temperature fluctuations were 
correlated (Equation (8)) as well as statistically inde- 
pendent (Equation (10)). They reported almost identical 
results. Hence, only the latter is shown here. 

The turbulent speCGtramQ@tor all thie’ runs were 
generally consistent and compared well with the -5/3 slope 
illustrated in Figure 25. There is c definite indication 
of a change over in the inertial subrange at a period of 
about 9 seconds (f = .11 Hz). This is to be expected 
since turbulence occurs at scales smaller than those 
associated with the waves from which turbulent energy is 
derived. in racao?eneaslope got thewspectra ase much nearer 
to -# than -5/3. Black (1965) has suggested this slopie is 
to be expected more often than -5/3 when there is density 
stratification. This would agree with the results sug- 
gested by the calculated structure function. These results 
appear to generally validate the assumption of statistical 


independence between turbulence and wave-induced 
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temperature fluctuations and that the natural preference 
in nature is for a near -5/3 dependence of turbulent 
spectral energy on frequency in the inertial subrange. 
Seitz (1971) has suggested that in a stationary, 
homogeneous, isotropic field of turbulence the energy 
spectra will have a dependence on a -2 power law at low 
frequencies and.a -5/3 power law at high frequencies pro- 
viding convection velocities are moderate. The transi- 
tion from one power law to the other should be at some 
intermediate frequency and would be a function of the 
magnitude of the convection velocity. Following this 
line of reasoning then it can be seen in Figure 25 that 
all runs exhibit inertial subrange character in the energy 
Spectra above 0.1 Hz. -Runs 1 and 8a, which had no 
apparent internal wave activity, exhibit similar slopes 
of —7/3 below 0.031 Hz with transition between these fre- 
quencies. However, Runs 6 and 7, which did show signifi- 
cant internal wave activity, exhibit slopes of" -3 and -4.3 
respectively below 0.031 Ha. Thus wouldgtuend Go reinforce 
the notion that internal waves contain much more energy 
Pome  lowersfrequenceresäand@are primariiy responsible for 


the time series being non-stationary. 








V. CONCLUSIONS 


The temperature field is extremely complex in the 
Shallow waters near the NUC Tower. There were indica- 
tions that double-diffusion type microstructure was 
present as well as the more normal, for this area, billow 
turbulence microstructure which appeared to be suggested 
by the temperature stratification for the majority of the 
runs. 

Temperature fields with no apparent internal waves 
may contain turbulent oscillations suggestive of internal 
waves or their remnants. The basis for this statement is 
thatthe osci lations in Run 1 appeared to We periodic Of 
period three minutes which is close to that reported by 
LaFond (1966) for internal waves. 

Significant internal waves present on the two days of 
the experiment were periodic with approximately 20 minute 
periods. In form they were more like solitary waves rather 
than sinusoids but this is probably due to the fact that 
they were shallow water waves at this location. Judging 
from the temperature profiles the internal waves appeared 
to have positive displacement from a mean position with 
little or no negative displacement. The positive displace- 
ments occupied less than one-half of one period. 

Internal waves appeared to be the major source of tur- 


bulence during this experiment under the conditions of 
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light winds and small ampitude waves. This conclusion 
is based on the higher variances in the runs with internal 
waves. 

When internal waves were present there was signifi- 
cantly more energy at all wavenumbers and not only at low 
vavenumbezs#, This statement iSi predicated on thel fact 
that even after filtering out waves of period 100 seconds 
and greater the runs with internal waves still had more 
energy than runs without internal waves. 

Internal waves may possibly break, thus contributing 
a large amount of turbulent energy to a localized area 
and thus assist mixing. 

With respect to spatial correlation lengths, the length 
of the array (64 inches = 1.6 meters) was too short to 
measure the scale length of the internal waves. The high- 
pass filtered time series had correlation lengths of 130 cm 
and less. There did not appear to be any depth dependence 
on the length but there may have been a dependence on 
whether the array was vertical or horizontal. 

With ragspect to the gradient field, it is confirmed 
that temperature gradients taken over short distances are 
exu emite 1) ata alb le Over *timerMand Psiplalce "as expected Mor a 
turbulenta phenomenon: 

Wavenumber spectra may be significantly altered at 
higher wavenumbers depending on the type of interpolation 
used in determining ‘the spatial correlation function to be 


HLourrer:  raoucrormed. Further, the vertical and perhaps 
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the horizontal particle displacements at the frequency 


of maximum coherence in the frequency spectra appear as 


changes of slope in the wavenumber spectra. 


The conclusions reached in respect -to the spectral 


medel proposed by MTMhorntonfet "al (1974) were? that: 


1. 


The model is not always realistic in the speci- 
fication of a constant mean gradient. The 
gradilente o requamed to besusicd “ameno ea ns Lex 
function to achieve consistency between the actual 
temperature spectrum and the wave-induced tempera- 
ture spectrum was much less than the observed or 
calculated mean gradients. 


The wave-induced temperature spectrum of the 
model showed great similarity to the actual tem- 
perature spectrum over a narrow band of frequen- 
cies associated with regions of high coherence 
betwen Surface waves and temperature spectra. 


The model describes reasonably well the wave- 
induced temperature fluctuations in shallow water 
(less than 18 m). Í 


The model illustrates that moderate waves do little 
mixing and serve only to pump the thermocline (and 
tenperaturei frend eu and donn: I = 


The temperature field lags the wave field by ap- 
proximately 180929; but the lag 415: modified by (as 
yet undetermined factors} and depth. ‘They spectral 
analysis yielded the above with respect to surface 
waves and the results of the isotherm plot in 

Run 3S Indicate the same for internal waves. 
TMeamouel demonstrates) the mhigh probability o£ 
Statistica 2ndependemees between turbulent and 
wass-inducedZtemperature fluctuations. 


The turbulent temperature spectra agrees very well 
xaıcıusınertral subrang theory in that the slopes of 
the 1og/log spyeetrazanie veuy near -5/3. That the 
slope is nearer to -2 may indicate that the time 
series was non-stationary. 


mnie lov frequency end of the turbulent spectra 
appears to have different slopes depending on 

whether the individual run had internal waves 

prasenits. 
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